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A computer simulation of a correlator receiver was 
developed and exercised to study the impact of a model-based 
signal processing algorithm on the detection of transmitted 
CW and LFM pulse acoustic signals incident on a planar array 
of electroacoustic transducers. The model of the ocean 
communication channel incorporates a space-variant sound 
speed profile. The transducer output electrical signals are 
cophased by an FFT beamformer via phase weighting, and 
summed to form a total array output signal. The total array 
output signal is correlated with a delayed replica of the 
transmit waveform and compared to a Neyman-Pearson thres- 
hold. Receiver performance is measured using a Monte Carlo 
technique to estimate the probability of detection for a 
fixed probability of false alarm versus the signal-to-noise 
ratio at the input of a single transducer. White, zero-mean, 
Gaussian transducer noise is assumed to facilitate compari- 
son between theoretical and simulated performance. Results 
indicate that model-based signal processing provides signi- 


ficant improvement of receiver performance. 
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Model-based signal processing is described by Mendel 
[Ref. 1] as an approach that exploits knowledge of the 
underlying physics of a problem to develop signal processing 
algorithms. Use of the approach implies that some a priori 
knowledge exists regarding the problem under consideration. 
In the case of an underwater acoustic communication problem, 
such a model has been developed by Ziomek [Ref. 2]. 

Ziomek derived a time-invariant, space-variant transfer 
function of the ocean volume based on the WKB approximation, 
which is an approximate solution of the linear, inhomogeneous, 
scalar wave equation describing the propagation of small- 
amplitude acoustic pressure waves when the speed of sound is 
a function of depth. Based on the transfer function of the 
ocean volume, Ziomek [Ref. 2:pp. 257-261] also derived an 
equation describing the output electrical signal at each 
element of a planar array of point sources. The output 
Signal is described in terms of the frequency spectrum of 
the transmitted electrical signal, the transmit and receive 
planar arrays, and the random ocean medium transfer function. 
Vos [Ref. 3] used these results to develop a computer program 
that generates time-samples of a real baseband output 
electrical signal at each element in a receive planar array 


as a function of a variable ocean medium sound speed profile, 


ıl 


planar array size, array far-field beam patterns and func- 
tional form of the transmit signal. Ziomek has since modi- 
fied this program to generate time-samples of the complex 
envelope of real bandpass output electrical signals. 

The research documented in this thesis has the following 
objectives: 

- develop a computer simulation of a correlator 
receiver which processes the output electrical signals 
generated by the computer program developed by Ziomek 
and Vos; 

- apply the concept of model-based signal processing 
to the development of the signal processing algorithm 
used by the receiver; 

- determine the effectiveness of the approach in the 
detection of signals from a planar array of point 
561176 0 ٦ 

Since the effects of the ocean medium on the signals 
processed by the receiver are embodied in the random ocean 
medium transfer function, the basic question to be addressed 
may be stated as follows. Can a priori knowledge of the 
ocean medium, based on physical principles of acoustic wave 
propagation, be used to improve the detection of signals in 
a receiver, and to what degree? 

Ziomek's use of linear systems theory to develop a 
transfer function model of the ocean medium immediately 
suggests the use of a compensating filter at the array 
Output to remove the undesirable time delays due to system 
geometry and wave propagation effects. This filter would 


ideally cophase the signals at each element in the planar 


array, resulting in maximum signal output when the signals 


1. 


are added together. This filter is implemented in the 
frequency domain through the use of Discrete Fourier 
Transforms (DFT), and the approach is exactly analogous to 
the FFT beamforming procedure discussed by Ziomek [Ref. 2: 
ES 153-1761]. 

In the frequency domain, the time delays due to system 
geometry and wave propagation effects are represented as 
phase shifts which may be eliminated if known a priori. The 
concept of model-based signal processing is applied here to 
obtain the proper compensating phase shift for the known 
system geometry and wave propagation conditions. 

Section II describes the theory used to develop the 
receiver model. The system context within which the 
SS ር Operates 1S described and related to previous 
investigations. A functional description of the receiver 
1s shown, and each of the major functional blocks is explained 
Bu some detail. Finally, a statistical description of the 
receiver's performance is developed. 

The computer implementation of the receiver structure is 
described in Section III. The logical flow of the computer 
program is discussed and related to the receiver descrip- 
tion. The use of multiple trials to estimate the probability 
of detection 1s explained. Each of the major subprograms 
is characterized in terms of function and implementation. 


Verification of the computer simulation is discussed last. 


ES 


Section IV presents the data obtained from the simula- 
tion when a rectangular-envelope, continuous wave (CW) 
pulse or a rectangular-envelope, linear-frequency-modulated, 
(LFM) carrier is transmitted. Receiver performance is 
described by plotting the probability of detection (Pd), 
for a given probability of false alarm (Pfa), as a function 
of the input signal-to-noise ratio (SNR) at each element 
in the receive array. Plots are provided for different values 
of Pfa, and show the relative improvement in receiver per- 
formance as various medium and wave propagation effects 
are compensated for by the model-based, signal processing 
algorithm. In each plot, the receiver performance predicted 
by theory when all array element output signals are precisely 
cophased, and the array element input noise is zero-mean, 
uncorrelated and Gaussian is shown as a dashed line. The 
dashed line is plotted from data obtained from a closed 
form expression relating Pd to Pfa as a function of array 
element input SNR, and is superimposed on the output plots 
to provide a baseline for judging the validity of the receiver 
simulation output data. 

Conclusions and recommendations are discussed in 


Section V. 
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POF THE RECEIVER MODEL‏ ط9۱ 


A receiver operates in the context of a total communi- 
cation system consisting of a signal source (transmitter), 
a signal propagation medium (channel) and a signal sink | 
(receiver). It is the model of the communication channel 
that is of initial interest since the signal processing 
algorithm will depend in large part on the physics des- 


cribing the propagation of the signal through the channel. 


OVERVIEW OF THE COMMUNICATION SYSTEM 

Ziomek's model (Ref. 2] of the ocean medium is described 
in general as a time-variant, space-variant, random filter 
(transfer function) in which the index of refraction, or 
equivalently, the speed of sound, is a function of depth, 
and includes both a deterministic and a random component. 
However, in describing the electrical output signals from 
the receive aperture, the model becomes more restrictive in 
the sense that the channel is considered to be time-invariant, 
but still space-variant. Furthermore, the transmit and 
receive apertures are taken to be rectangular, planar arrays 
whose elements consist of complex weighted point sources. 
Complex weighting of the array elements provides the means 
for amplitude shading and beam steering both the transmit 


and receive array patterns. The complex weights are ideal 


o 





for removing the undesired effects of the channel engine 
output electrical Signals inom the receive array Ten ne 
and become the tool for applying the model-based, signal 
processing concept. Figure 1 depicts the geometry Gf tne 


transmit and receive arrays in the transmission medium. 
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Figure l. System Geometry 


Ziomek goes further than simply considering system 
geometry in generating the output signals at the receive 
array. The spatial variance of the transfer function per- 
mits modeling of propagation effects due to the dependence 
of the index of refraction on depth. Figure 2 shows the 


ray-bending effects of propagation through such an 


l6 
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Figure 2. Ray Path Bending Due to Inhomogeneous Medium 
inhomogeneous medium. The phase shifts due to system geometry 


and propagation effects may be considered together or 
Separately in the generation of the output electrical signal 
data in the computer simulation due to Vos [Ref. 3]. It 

is this output electrical signal data which is used as the 


input signal to the receiver. 


Emer UNCTIONAL DESCRIPTION OF THE RECEIVER 

The output electrical Signal data processed by the 
receiver are time samples of the complex envelope of the 
bandpass electrical signal at the output of each element in 
the recelve array. This implies that the bandpass acoustic 


Signal incident on the receive array has been converted to 
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an output electrical signal by each transducer element in 
the array. Each real, bandpass electrical signal then 
passes through a ә demodulator to become an 
equivalent baseband, complex envelope signal, which is 
time-sampled and converted from analog to digital form. The 
complex envelope is represented by the I-channel (in-phase) 
and Q-channel (quadrature-phase) components generated by the 
quadrature demodulator.  Time-sampling is done in a manner 
that satisfies the Nyquist criterion for the baseband infor- 
mation contained in the I and Q channels. Thus, many of 
the components associated with a receiving system are already 
contained within the simulation that generates the planar 
array output signal data. The receiver simulation assumes 
these components exist, and essentially provides signal 
processing of the complex envelope of the output bandpass 
electrical signals. The major functional blocks of the 
receiver include: 

- an array signal processor, 


- a correlator implementation of a matched filter 
receiver, 


- a magnitude-square operation, and 

- a threshold decision operation. 
Figure 3 shows the functional block diagram for the receiver 
model. The magnitude-square and the threshold decision 
Operations are considered to be part of the correlator 


(matched filter) detector function block shown in Figure 3. 
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Етик 5. 


Receiver Block Diagram 


The receive array is assumed to be a rectangular planar 
array of MxN elements where both M and N are odd numbers. 


Each element is assumed to be an omnidirectional point 


፡ከ ለ ን. The geometry of the planar array and associated 


mathematical notation is shown in Figure 4. 
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Figure 4. 


Planar Array Geometry 
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The function of ecach'"usadraturesdemoduldtore ән 
convert the amplitude and angle-modulated, bandpass, 
electrical siqnal at cacızeraneducer “06007 Ди ms 
baseband, complex envelope. Thus, the low-pass complex 
envelope may be sampled at a much lower rate. Note that the 
use of a quadrature demodulator does imply that the carrier 
frequency of the transmitted signal is known. It should be 
emphasized that the output electrical signal at each element 
in the receive array is passed through its own quadrature 
demodulator before array processing begins. The quadr omi 
demodulator is shown schematically in Figure 5, where the 


complex envelope of y(t), denoted y(t), ys given by, 


y(t) = y (t) m ነ (25 
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Figure 5. Array Element Quadrature Demodulator 
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where: 


y 1777177 75 UI Channel" in phase) component 
6 
2557 ana 
E is the O-channel (quadrature-phase) 


Component Of y(t). 


The array signal processor, shown schematically in 
Eure 6, 1S a FFT beamformer. The function of the array 
processor 1S to maximize the total output signal when the 


Signals from each element in the array are summed. 


(ITs, md, dy)  Y(9f md,,nd,) 
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eure GO. Array Processor Block Diagram 


The array processor is essentially a filter that compen- 
sates for system geometry and wave propagation effects on 
the signal transmitted through the channel. The array 


processor is implemented in the frequency domain where 
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filtering can be obtained by multiplication of each spectral 
component of the complex envelope by a complex weighting 
coefficient. The complex coefficients are computed from 
a knowledge of the channel model, thus the application of 
model-based signal processing. After applying the complex 
weights, the Inverse DFT is computed to recover cophased, 
time domain, complex envelope signals from each array ele- 
ment. By summing these cophased signals over all M*xN array 
elements, a constructive interference effect is achieved, 
and the total output signal is maximized. 

The matched filter portion of the receiver is imple- 
mented by correlating the total time-sampled output signal 


and noise from the array, A Where, 


፲ (5፻7) = yy (2TÇ) + ng GT) ( 2:599) 


with a time and frequency shifted replica of the complex 
envelope of the transmitted waveform, x(t). Since the 
phase of the received signal is, in general, unknown, the 
magnitude-square of the correlator output is taken as the 
input to the threshold detector. This input is compared to 
a preset threshold level y to determine the presence of a 
signal. The preset threshold y is computed from a Neyman- 
Pearson criterion. The schematic of the correlator/matched 
filter detector is shown in Figure 7, and has been shown by 


Van Trees [Ref. 4:pp. 244-247] to be the optimum receiver 
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۳-7607. Correlator/Matched Filter Detector 


for the detection of a bandpass signal, with random ampli- 


tude and phase, in the presence of white, Gaussian noise. 


መጠበ  IISTICAL DESCRIPTION OF THE RECEIVER 

The output electrical signal data embodies several 
assumptions and restrictions that must be restated before 
further development of the receiver model. The signal data 
is assumed to be generated by a real, bandpass, acoustic 
field, y,(t,r), incident upon an array of ٠166 ETC 


transducers. The acoustic wave, ۴۴٢ح 9و‎ 


YM 


^^ 


in the የ“ directson vith ህፎ 1ርነፎሮ በ6962. or 
yuy(t,r) - x(tilr:njl/c) 28) 


where; 
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г = xx + yy + Zz )2 


Mo PS + ex. + w. (2250 

sin 9 cos Yo (2.6)‏ کت 

və s sin 00 Sani Vo Фи 
апа 

wy = Cos 8j (2.8) 


Note that x(t) may be an arbitrary function of time. 
l. Array Element Output Signal Description 
The acoustic field, incident upon an element in the 
array, is converted to an output electrical signal by the 
transducer, and is transformed into a baseband complex 
envelope at the output of the quadrature demodulator. After 
time sampling at a rate f = 1/T. (in samples per second), 
Spatial sampling over the receive planar aperture in incre- 
ments of a, and .. and assuming linearity in the transducer 
operation, the complex envelope of the output electrical 
th 


signal at them ined element of the planar array lying in 


the XY-plane may be written as, 
" ۵ = X (£T. y [ugmd., * .... (250 
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For example, if the real, bandpass output electrical 


10) 


Sn 


17 


21) 


.14( 


15) 


signal y(t) is an amplitude and angle modulated cosine wave 


ا 


the complex envelope of y(t) 


—. 


of the form 


y(t) - a(t)cos[2Tf t 5771 


A ۶۱ 


al cos 00) 


a(t)sin 6(t) 


/ 2 
y^ (t) + لا‎ 


rc) 
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after quadrature demodulation, 


may be represented as 


where: 

yg, (t) 
and 

yore) 


or, in magnitude-phase form 


alt) 


and 
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It is also possible to show that the envelope function, 


E(t), of the real, bandpass signal y(t) is 
E(t) = la(t)l (2265 


The main advantage of vorking vith the complex 
envelope form of the real bandpass signal is that the 
baseband complex envelope may be sampled at a rate Eo (in 
samples per second) determined by the bandwidth of the 
baseband modulating waveform, independent of the carrier 
frequency. 

It is further assumed that the complex envelope of 
the transmitted waveform m can be represented exactly, 


over an interval of 1 seconds, by a finite, complex, 


Fourier series such as, 


~ K 
x(t) - 2. c,exp[*j2maf,t] (220666. 
qe-K 4 


vhere the complex Fourier coefficient 8 can be written as, 


с = +36 2 
с 7:7: ) q! ( ) 


where: 


f_: 15 the fundamental frequency in Hertz of 
the srgnal x(t) with period + 

0 0 
seconds, and 
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K: is the maximum number of harmonics used to 
represent the signal x(t). 


The complex Fourier coefficients Ca can then be 


determined directly from the DFT with respect to time of 


2... o 
I 
1 = eue 
c = = } x(2T_)W (2.19) 
q L@ L=-L' S L 
where: 
A 1 a 
Wr =  exp[-*j2m/L] 1220) 
‹ ፏ 
L = (va bee (55 
and: 
p is the total number of time samples taken 


during the time interval T, seconds where 
L is a non-negative, odd integer, 


Tj: is the fundamental period or data record 
length in seconds, and 
T. : is the sampling period in seconds (note 
S 


TO satisfy the Nyquist sampling theorem, it can be 


ol [Ref. 2:pp. 164-165] that, 


)5 710 و7 
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and, 


Poes mo (2: 2 


2. Generation of Complex Weight Phase Factors 
With these assumptions it can be shown [Ref. 2:pp. 
165-166] that the normalized DFT with respect to time of the 


complex envelope of the output electrical signal لا‎ )1 md r 


nd.) can be expressed as, 


~ 


25750 - cg exp [*39,]expIzj2nfu m, /c] 
'"* exp [336, J exp Ix 32mfv nd /el (2.22) 
where: 
f = a + qf. 5/50 E (25 
ana: 
Ед: is the carrier frequency, 
C : 15 the constant speed of sound in homogeneous 


Ocean medium at the receive array, 


E is the direction cosine of the wave propa- 
gation vector along the x-axis, equation 
(2:6 ame 

v : 15 the direction cosine of the wave propaga- 

0 ; ን : 
tion vector along the y-axis. See equation 
150 
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Phase factors ди апа ቀ; are due to the separable 
complex weight, Cnn" That is, 


C - с̧а E Am Xp 1+30, 1b,expl+309,1 (2.26) 


ከ ከከ that by a proper choice of Ön and фи in equation 


(2.24) where, 

0 = + 271211107 © (257) 
and 

ф = ስ (2220) 


the phase shifts due to system geometry may be completely 
cancelled leaving only the complex Fourier coefficients c 
at each element in the planar array. It should be noted 
that the phase correction factors are functions of both the 
array element position (mdx or ndy), and the frequency f of 
each spectral component of the input bandpass signal spec- 
trum. That is, the frequency is given by equation (2.25). 
Once the phase shifts due to geometry are eliminated, 
taking the inverse DFT with respect to frequency of the 
Spectrum of the output electrical signal Y (gf, md, , nd, ) 
will yield cophased time signals at each element.  Summing 


all the cophased signals will result in a maximum total 
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output signal Yip (LT) from the array processor. That 


15 , 
Ë M' ፐ| ር 
ای‎ acr a. y (AT, md, ,nd.) (22200 
where: 
M" = (M-1)/2 (2.30) 
and 
N' = (N-1)72 (2998 


This approach is extended to compensate for the 
deterministic signal phase shifts caused by transmission 
through an inhomogeneous ocean medium. Since the variation 
in the speed of sound c(y) is assumed to be a function of y 
(depth) alone, the calculation involves only the Que phase 
factor of the complex weight. Ziomek! has shown that for a 
sound-speed profile c(y) with a constant gradient g, a 
closed form expression for the deterministic component of 
phase shift is, 

MC 


E O 
Фмр pO ٠٦ ےت‎ np?) -1] + Ay (245 


Extension of work by Ziomek based on expression for 9 MD 
[Ref. 2:pp. 263-268]. 
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5) 


. 34) 


75: 


550) 


np(y) = co/epiy) (2 
Cp (Y) کو تہ‎ улу د‎ (2 
Лу = ви اج‎ Us me (2 
Ky = 2170 (2 


is the space-variant (with depth only) 
index of refraction, 


is the speed of sound at depth y, 


is the speed of sound at the transmit 
array, 


is the wave propagation constant at 
the transmit array, 


is the gradient (slope) of the sound 
speed profile, 


is the depth of the center element of 
the transmit array, and 


is the depth of the center element of the 
receive array. 


where: 


Yr 


The negative of the deterministic medium phase 


(2.32), is simply added to the system 


0-202, equation 


geometry phase factor, equation (2.28), to obtain the total 


complex weight phase factor in the y direction. 
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3. Array ProOcessorWOUtonENOTSghU NEST SED SS 

Development of a receiver model also requires a 
Specification of the noise environment in which the receiver 
Operates. Any realistic model for the noise environment is 
extremely complex. However in order to form tractable 
theoretical results that can be reasonably approximated in 
a computer simulation, zero-mean, additive, white, Gaussian 
noise (AWGN) is assumed at the output of each array element. 
Remember that the array element output is also the input to 
the quadrature demodulator. 

The AWGN model permits derivation of a closed form 
expression relating Pd, Pfa and array input Signal-to- 
Noise Ratio (SNR). This type of noise process can also be 
reasonably approximated by uSing computer generated pseudo- 
random number ም ans from a standard Gaussian random 
number generator. By comparing the theoretically predicted 
receiver performance with the results of the simulation, 
verification of the computer implementation of the receiver 
can be achieved. Once verified, the computer simulation 
can test more realistic noise models with some confidence 
in the resulting data. 

The input SNR at a single element in the array is 


defined as, 


a 2 
E{ ly (4T, ,m. , nd, ) | } 


(сә 


SNR; (2300) 


E 2 
E( |n (4T, , mà. , nd. ) | ) 
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where Ei } denotes the expected value of the quantity within 


LS ergodic, 
y! - 


the curly braces. 


If the random input signal y (4T, , md, , nd 


then the mean-square value of the signal can be found by 


computing the time-average instead of the ensemble average,. 


B 2 
«|y (AT, md. nd, | » 55: 


time-average of the quantity vithin 


T 2 


dt (250) 


equation (2.39) can be approximated 


purposes by, 


(2.40) 


—u (2 11) 
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that is, 


Y 2 
EC|y (27, md. , nd.) | — 
where -(:)» denotes the 
the parenthesis, or 

«(2 x 
0 ۔‎ 72 


The integral in 


for computer simulation 


mere the infinitesimal 


Because of the assumption regarding the use of a 
finite Fourier series to represent the original transmitted 
waveform, the input signal power at a single element in the 
receive array can be computed from the sum of the magnitude- 
Square of the complex Fourier coefficients [Ref. 5:pp. 44- 


857 ۲135۲15 


бу 2 К 4 
“Ју (Отс̧ула,па,) | . ) eq (22፡8 
q--K 
These coefficients are easily obtained for each array 
element by computing the DFT with respect to time of the 
complex envelope of the output electrical signal data at 
each element, or 
Ies ~ 
EN XL 42 
7 = 2 y (QT mà. ,nd )Wr (2.43) 


The mean-square value of the zero mean noise signal 
at the array element input is equal to the variance a or 


the noise input, or 


~ 2 2 
E( [n (2T, r ) | } = o (2.44) 


Substituting equations (2.42) and (2.44) into equation (2.37), 


the expression for input SNR at element (m,n) may be written 
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35 , 
SNR, = NN CIIM D (2545) 
or, by rearranging the variance, equation (2.45) becomes 
O = E |^ /SNR (2.46) 
و‎ 1 l 


The measured input signal power at an array element 
and a desired input SNR parameter value SNR; can be used to 
obtain the noise power (variance) required to scale the 
output from the random number generator. The ability to set 
a desired input SNR value is necessary in order ts test 
receiver performance. 

Because the DFT and IDFT are linear operations, the 
noise statistics at the output of the array processor are 
Still Gaussian, and uncorrelated in both spatial and 


temporal coordinates. The total noise signal ng CAT.) at 


the output of the array processor may be written as, 


: M' N' 
MED - ) )  n(£T,,md, ,nàd, ) (2.47) 
= . m--M' n--N Ы = Y 


; | ду. 
The variance of the total noise Om 1S equal to the mean 


Square value of the total noise signal from the array 
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processor, and is also equal to the sum of the variances 
of the noise signals at each element since the noise process 
at one element is by assumption independent of the noise 
process at any other element in the array. The variance 


of the total noise signal is then 


(In 2) M' N' 2 
E (LRI) = O (2.48) 
T = | i acia = 


: 2 
I£ the variance DE at each element is the same for all 


elements in the array, or 

O = О (2249) 
then, the variance of the total noise o: can be written ae 
terms of the variance of the noise at each element in the 
array by substituting equation (2.49) into equation (2.48) 
EO Əbül. 

O —  MNO (2095099 


Or, in terms of mean square values, 


E{|n,(2T,)|°} = MNE{|n(27,) |*} )2.51( 


DIG 


where E{|n(2T_) |*} is the noise power at an array element 
when the noise power is considered to be the same for all 
array elements (m,n). 

Both the output signal and output noise components 
at all elements in the array are baseband complex envelope 
signals with I and Q channel components. The bandwidth of 
the noise signal is set by the bandwidth of the low-pass 
filter in the quadrature demodulator to W Hertz. For the 
purpose of the simulation, W is always adjusted to include 
the highest harmonic of the complex envelope signal, that 


is, 


W = Kf (2552) 


Knowing the bandwidth W and using the AWGN assumptions, 
the noise power spectral density be at each element can be 
related to the variance of the noise process at each element 


By, 


N = : 


- mn (2K£,) (255) 


The SNR at the output of the array processor SNR, 15 


defined as the ratio of the total signal power to the total 
noise power, or 
E(1y, (672) 17) 
A T S 


SNR fO 6. (2.54) 
arb ٠٦ 
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where, 


95 2 ۵ l 
ELİy, (Um) ис ии) 1 y (4T. md. , nd.) 
m--M' n-z-N' 


[n (2.55) 


For the case of perfectly cophased signals with identical 
signal power at each element, the right-hand side of equa- 
tion (2.55) 67 0 


2 


EC] y, GT 17) = (MN) EL |F OTL) | ) (2.56) 


where EC|y (17, |^) is equivalent to the time-average power 

at each element in the array, and 8 ٥٣٥٦ 15 equivalent 
to the time-average power of the total output signal from 

the array processor. Using the uncorrelated and equal 

array element input power noise assumptions, the expression 
for SNR, can be rewritten by substituting equations (2.51) 


A 
and (2.56) into equations (275 Es 


5(|ሃ(፤፲.) |") 
SNR, = (MN ———— = MNSNR, - (2.57) 
Ej ipm P 


Since the array gain (AG) is defined as, 


AG 10 log, , ISNR,/SNR, ] (dB) (2.518) 
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the array gain in this case becomes, 


AG = 10 log, q IMN] (dB) (2.59) 


eé that fOr a 5x5 element planar array, M*N = 25, and 
the array gain is 13.979 dB. 
ENNI LSbhesrs restsng and the Neyman-Pearson Criterion 
The correlator/matched filter detector portion of 
the receiver, Figure 7, is modeled as a binary hypothesis 
testing problem using the Neyman-Pearson decision criterion. 


are defined as, 


The two hypotheses, Hy and Ну, 
r (2T,) = (25 ии 
nn (TÖ : E e 


where, in terms of the transmitted waveform x(t), 


M' N' . 
lI) = c uM is x [ (4T Tmn? ema, nd] 


x exp [*5270, 7T .] (25: 
and ит) is given by equation (2.47). If one assumes 
that all array element output signals are cophased and 


identical, equation (2.61) reduces to 
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~ 


Yo GT.) = CMNx (L£T_-t,) exp[+j27o,2T.] (29995 


where: 
C = a exp [+30] (2462, 


апа, 


T.: is the actual time delay in seconds at 
each element in the array due to range 
separation between individual elements 
of the transmit and receive arrays, 


is the actual time delay in seconds due 

to range separation between transmit 

and receive arrays when all signals are 

cophased, 

is the actual: doppler shift in Hertz due 
to relative motion between transmit and 

receive arrays, 


a: is the amplitude attenuation factor that 
is, in general, a random variable, 


5 is the generalized phase shift of the 
received signal with respect to the trans- 
mitted signal. (In general, it is also a 
random variable dependent on both spatial 
and temporal coordinates.), and 

Ho: is the null or noise only (no signal) 
hypothesis. 

Matched filtering is obtained by correlating the 
complex envelope of the total array output r(2T_) with the 


complex conjugate of the processing waveform UN The 


functional form orf g(2T,) can be written as, 


g(2T_) = (QT -1) exp[+52ng2T 1 (2.64) 
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which is a time and frequency shifted replica of the 
transmitted signal x(t) where 
т: is the estimate of time delay at the receiver, 
and 
$: 15 the estimate of the doppler shift at the 
receiver. 
For the purposes of this study, the estimates of 
time delay and doppler shift are assumed to be precisely 


5፡8 22, 0 


H 5 6) - 0-0 (2.65) 


and 


x00 (2:66)‏ ت2 ФА‏ م 


መሪ 


The correlation, or inner product £ between two 


munetions r (t) and g(t), is given by 


+00 
Ie (2.67) 


I> 


2 = <r(t),g(t)> 


which is approximated in the simulation using the trapezoidal 


male approximation to the integral, that is, 


| L'-1 ሯ (5፲..) 9* (5፲. ] + 2[(፤+3)፲ 19* [(4+1)፻. 
405 ا‎ Ми —— = с (2:55: 
=-1, S 
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The magnitude-square of the correlator output is 
taken for two reasons. First, the phase of the carrier 
frequency waveform is, in general, unknown. As the phase 
of the received carrier varies with respect to the phase of 
the quadrature demodulator local oscillator (LO) signal, the 
output of the quadrature demodulator would vary from a 
maximum negative to a maximum positive value depending on 
the phase difference between the carrier and the LO which 
usually is taken to be a uniformly distributed random varia- 
ble between 0 and 27 radians. The change in polarity of 
the quadrature demodulator output would propagate through 
to the output of the integrator in the correlator/matched 
filter detector, Figure 7. Taking the magnitude of the 
integrator output ensures that the input to the threshold 
comparator will always be non-negative regardless of the 
phase difference between the carrier and LO waveforms. 
Second, when the array element noise statistics are assumed 
to be Gaussian, the square of the magnitude of the integrator 
output yields an input to the threshold comparator which 
can be described statistically by exponential density 


functions for both Hy and H, signal hypotheses. As will 


0 
be shown in the following derivations, the exponential den- 
sity functions can be easily integrated to obtain a closed 
form expression for Pd and Pfa in terms of the SNR at the 


input to the threshold comparator. The output of the 


magnitude-square operation is the sufficient statistic 
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E hich the binary decision is made. That is, choose Hy 


ME, 


مد 


5 (2.69) 


or choose Ho If, 


2 


a 75750) 
اا‎ ርን 


Assuming that the total noise n, (t) is a baseband 


1 
Gaussian process of bandwidth W, the conditional probability 
density functions (pdf's) of the magnitude-square of the 
correlator output with and without a signal present can be 


shown to be 5: The conditional pdf"s are given 


by 


EN 1/ (205) exp L- | £| ^7 (2051 (225710 


and, 


р(|2|“)н|) 1/ (20%) exp[-|2£|*/ (204) ] 10002) 


For zero mean, AVGN, and cophased, equal eneroy (pover) 


signals at each array element output, the variances of 


“Derivations for the pdf"s, expressions for Pd, Pfa and 
the decision threshold Y were provided by Prof. L. J. Ziomek 
iu private communication. 
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the magnitude-square of the correlator output can be shown 


to be 
o? 592772 5 P (2.73) 
O N 0 A | 
and 
2 2.17 22 
oq “ELİ ası. 
2 2 2 
= (MN) E )|Х(т/9) | + 0ط‎ 32 (2.74) 
where: 
SR 

EEN : is the magnitude-square of the correla- 
tor output when the input to the 
receive array is noise alone, 

|£ |^: is the magnitude square of the correla- 
tor output when the input to receive 
array consists of signal alone, 

N is the power spectral density level of 

O 
the noise signal at each array element 
ت6‎ lir: 

F. is the energy in the transmitted signal 
which for simulation purposes is defined 
to be equal to the energy Ez in the local 
processing waveform g(t), 9 

2 . 

Ela ]: is the mean-square value of the amplitude 
attentuation factor 3 7 i 
only when deterministic effects are 
considered), and 

2 . ፡ 
IX( 1,0) | : is the magnitude-square of the auto- 


ambiguity function. Note that İX((r,9) 15 = 
Е2|Х,(т,ф) |2, or in our case where 1 - 0 
and ¢ = 0, |x(0,0)|2 = EZ since 

- 1. ГВеТ ил ерр 707 


о 
о 
NO 
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TO calculate the threhold y for a desired Pfa and 
12 oNR, the integral of the pdf for the null, or Ho 
hypothesis, is set equal to the desired Pfa, and the result- 
ing equation is solved for the threshold. The magnitude- 
square operation results in particularly simple threshold 
equation when the input noise is assumed to be white, zero 


mean, and Gaussian. Computing the Pfa yields 


+00 Es 
ə 2 
/ (iz. JH) dle | 


Pfa 
1 


exp [-v/ (205) ] 


exp[-y/(MNN Es) ] (2275) 


By solving equation (2.75) for the threshold y, 
Y - MNN,E. in [1/Pfa] (2.76) 


Once the threshold is obtained, the correlator 
output pdf for the H; hypothesis may be evaluated to obtain 


the Pd. 


OO 


nr 2 
7: Hy) ale) 


Fd 
1. 


exp [-y/20%] 


exp(-y/1 (MN) “E(aYİx (r, 6) [I RMNN EZ) (2.77) 


45 


Since in our problem Tt = 0 and 6 = 0) 


2 2 2 
IX (1,9) | = İX(0,0)) E (2:72 
and if ve let 
2 
Ela JES 
ик Ne E 2.79 
Wa N- ) ( ) 
then, 
2 
MNE{a JE; 
Pd = expi-y/ [MNN, E; ( N + 1)]) (2.80) 
anda 
Pd = expí-y/[MNN E-(SNR,.,2 + 1)]) (22588) 
DUZ e F 


Substitution of equation (2.75) znto equation (299 ٤٣٦ 
gives the desired closed form expression relating Pd, Pfa 
and the SNR of the magnitude-square of the correlator output, 


that is, 


Pd = Pfa (2.92 
This result agrees in form with the result given in Van 
Trees [Ref. 4:pp. 246-247] for a similar single channel 
receiver model. Figure 8 graphically depicts the relationship 
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Figure 8. Density Functions of the Magnitude- 
Saler cor relator Output 


between the conditional pdf's, the decision threshold y, the 
Aia and the Pd. 

iln corroluator olüEput SNR, equatıon (2. /9), can be 
related to the input SNR at a single element in the array 
through the array gain, and a factor resulting from the 
255: ly different definitions of input SNR and correlator 
output SNR. In terms of array element input signal energy 
Es; the input SNR given in equation (2.45) may be rewritten 


as, 


(1/T ¿ES E(a” )B- 
15 WA :  ... (2.83) 
G TAO 
mn o mn 
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where, 


_ 2:5 
ES -0 a JE (2.22) 


is the average received energy at a single element in the 
array due to the transmitted signal x(t), and from equation 


155) 





0 = 2kf N, = 2 (2.85) 


The magnitude-square correlator output SNR can be written 


as, 
E{|2,|°} 
F 
where, 
E{|2.|°} = (MN)2E(a JE (28፡91) 
and, 
EC|Z Ie = MNN_E~ (2.88) 
N O x 


Substitution of equations (2.87) and (2.88) into equation 


(2.86) gives 
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"7 
— X 


N 
0 


SNR 2 = MN 2) 
iE | 
Pyesubstituting equation (2.85) into equation (2.83) and 
rearranging terms, equation (2.83) becomes 
2 
LX )E- 
SNR, = اق‎ ሀ” (2500 


ZKN 
0 


By rearranging terms in equation (2.90) and substituting 
the result into equation (2.89), the desired relationship 
expressing the magnitude-square correlator output SNR in 
terms of the array element input SNR at a single element 


in the array is obtained. That is, 


700 ۴ = MNZKSNR; (2. Ли 


Writing both sides of equation (2.91) in dB form yields, 


|| (dB) = 10 log 109 (MN) + 10 log (2K) + SNR; (dB) 
“an” 


AG 12502) 


where AG is the array gain given in equation (2.59). 

To summarize, Section II.C provides the equations 
needed to implement a computer simulation of the receiver 
structure described in Section II.B. The implementation of 


the computer simulation is the subject of Section III. 
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Tit. COMPUTER Ss 1 سن‎ OE TNE TECET EE 


The computer program RCVR simulates receiver operation 
through straightforward application of the equations and the 
concepts developed in Section II. Written in FORTRAN, the 
computer program RCVR consists of a top level, or main 
program, and nine subprograms. The computer program will 
be explained from a functional viewpoint. That is, the 
algorithms used to implement the receiver simulation will be 
related to the theoretical development outlined in Section 


II, but translation of these algorithms into FORTRAN state- 


ments will not be discussed. The main program will be 
described first. The description of the main program will 
be followed by a detailed discussion of each subprogram. In 


addition to explaining the computer simulation, the methods 
used to validate the receiver simulation output data will 


be presented. 


A. TOP LEVEL PROGRAM DESCRIPTION 
The organization and logic flow of the top level, or 
main program, is shown in Figures 9a through 9e. The 
functions of the main program include: 
- initializing the simulation run-time environment, 


- invoking subprograms in the proper sequence to 
process the input signal data. 


- providing control logic and noise generation algorithms 
needed to measure receiver performance, and 
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Figure 9a. Program RCVR Flowchart 
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Figure 9e. Program RCVR Flowchart 


- selecting the form of the simulation output data. 
Initialization of the run-time environment involves reading 
input data from file to internal storage, generating the 
local processing waveform, measuring the time-average input 
signal power at each array element, computing the baseband 
signal bandwidth, and generating the complex weights used 
by the array processor. The input data is read by a call 
to the subprogram READY which returns a set of simulation 
parameters in COMMON storage and the complex envelope 
electrical signal data for each array element in matrix 
form. The local processing waveform is obtained by a call 
to the subprogram SGNLGN. The time-average signal power 
at each element (m,n) in the planar array is found by using 
the fact that the original transmit signal vas SNE | 
from a finite Fourier series, and applying equation (2.43) 
to obtain the complex Fourier coefficients, Ван Once the 
complex Fourier coefficients at each element in the planar 
array are obtained, the time-average power is computed using 
equation (2.42). The baseband signal bandwidth W is com- 
puted using equation (2.52). The separable complex weights, 
Cm and d, inequation (2.26) are generated in two stages. 
First, the real amplitude factors, ag and b, in equation 
(2.26), are computed by a call to subprogram AMPWGT. Seconde 
the real phase factors, om and on in equation (2.26), are 
obtained by a call to subprogram PHSWCT. The complex weights 
are then computed by combining the amplitude and phase 


factors as shown in equation (2.26). 
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The first signal processing step is to generate the 
total array output signal in the absence of additive noise 
07 a Call to Subprogram ARYPRO. Subprogram ARYPRO which 


uses the complex weights c 


m and др to cophase the planar 


array output signal returns a total array output signals. 
The total array output signal energy is computed using 
equation (2.40) for later use in calculating a SNR at the 
output of the array processor. The ratio of array output 
SNR to the input SNR defines the array gain, and provides a 
check on the validity of the data generated by the array 
processor algorithm. 

Receiver performance is measured by computing a relative 
frequency estimate of the Pd over the specified range of 
array element input SNR values when the Pfa is a known, 
time-invariant parameter. Using the Pfa and a range of SNR 
values specified by the programmer in a receiver control 
data block, the main program establishes nine input SNR 
values for which the estimate of the Pd will be computed, 
and determines the number of trials, or runs, to be used in 
computing the estimate. 

7 Point E in the flovchart of Figure 9c, the main 
program enters a loop that initializes the noise source, 
performs an array “kini calculation, computes the detection 
threshold and then enters a second, inner loop vhere the 
correlation detection is done. After exiting the inner 


loop, the relative frequency estimate of the Pd is computed. 
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The outer loop increments through the input SNR values 
determined earlier, and exits the loop when the specified 
maximum input SNR value is reached. 

The first step in initializing the noise generator is 
to determine the noise variance (power) at each array element. 
The power level, or variance, of the noise samples for a 
particular element in the array is a function of the time- 
average signal power at each element and the array element 
input SNR. The variance can be found using equation (2.46). 
By controlling the power level of the noise process at 
each array element, the receiver performance can be measured 
over a range of array element input SNR values regardless 
of the time-average signal power level at a particular 
array element. 

The variance computed in this manner is the variance of 
the complex envelope baseband noise signal. The variance of 
the real baseband noise signal produced by the noise generator 
subprogram AWGN in the I or Q channel is one-half the vari- 
ance of the complex noise signal since the complex envelope 
is the sum of the two independent I and Q channel noise sig- 
nals. Therefore, the variance of the complex envelope of 
the noise is divided in half prior to being passed to the 
noise generating subprogram as a scale parameter. 

The variance of the total complex envelope noise signal 
oz is obtained by summing all the noise variances of the 


“b 


! 2 
array element complex envelope signals, o 


mn” 88 indicated 
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by equation (2.48). A power spectral density of the total 
noise signal is computed using equation (2.53) where oA is 


substituted for ed 

The SNR at the array processor output is computed using 
equation (2.54) where the time-average power of the total 
array output signal is substituted for the mean square 
ensemble average, and the mean square ensemble average of 
the total noise is taken to be the same as the total noise 
signal variance. The array gain is computed using equation 
(2.58), and is held in internal storage for later output in 
EB ular format. 

The detection threshold y is computed using equation 
(2.76) vhere the total noise pover spectral density is 
substituted for the MNN. factor and the energy of the local 
processing waveform E is used instead of the transmit 
Signal energy E-. For the purpose of this simulation however, 
the energy in the local processing waveform is the same as 
the energy in the transmit waveform. That is, 52 and ES 
are equal. 

The second, or inner loop, of the main program begins 
E oulnt D in Figure 9c. Within the inner loop, the complex 
envelope noise signals are generated for all array elements, 
the total noise output signal is generated by the array 
processor subprogram ARYPRO, the total signal and total 
noise are summed and correlated with the local processing 


waveform, and the signal detection decision is made. The 
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inner loop terminates when the number of iterations through 
the loop equals the number of trials allowed to form the 
estimate of Pd. 

The time-sampled, complex envelope array element noise 
data is obtained by repeated calls to the noise generator 
subprogram AWGN. Each call to AWGN returns a properly 
scaled pseudorandom number representing one sample of the 
noise process in the I or Q channel at a particular array 
element with L time samples taken for each of the I and Q 
channels. The noise signal data is stored in the same 
matrix form as the array element output electrical signal 
cara 

The complex envelope noise data is submitted to the 
array processor subprogram for processing in exactly the 
same manner as the input signal data. The result is an 
array total noise output signal. Processing the noise alone 
in this manner provides some gain in execution speed and 
provides the flexibility to estamacequbespobsl qug ٣ 
false alarm directly, if desired. The DFT and IDFT are 
linear operations, and the principle of superposition holds. 
Therefore, the addition of the total noise and total signal 
at the output of the array processor is equivalent to adding 
noise to the signal at each element prior to passing the 
data to the array processor subprogram. 

Correlation of the sum of the total signal and total 


noise with the local processing waveform is accomplished by 
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using the trapezoidal rule to approximate the integral of 
equation (2.67). The magnitude-square of the correlator 


output is obtained by taking the complex product 
- 3 it 
| 2| ис, )0 13] 


The threshold detection portion of the receiver is 
implemented by comparing the output of the magnitude-square 
operation to the decision threshold, using a simple IF-THEN- 
ELSE binary branch. The number of hits, or times the output 
exceeds the threshold, are counted and stored. The process 
of generating noise samples, and making a hit or miss deci- 
sion continues through a large number of trials. Since the 


correlation is done using the H, (signal plus noise) hypothe- 


1 


sis, Pd can be directly estimated using the ratio 
Pd = HITS/TRIALS ጋ) 


The first approach taken to determine the minimum number 
of trials required to estimate Pd was based on equation 
(2.82). After rearranging terms, equation (2.82) can be 
written as 


pg | 1 * SNR] 


Pfa 5 


The central idea vas to compute an estimate of Pfa, using 


equation (3.3), from the relative frequency estimate of Pd 
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in equation (3.2). The algorithm would then terminate when 
the computed estimate of Pfa differed from the Pfa input 
parameter by some arbitrary small amount. However, the use 
of equation (3.3) was found to be a poor test for establish- 
ing when the algorithm should terminate, and would not 
terminate the algorithm for Pd values greater than about 
0750: 

A perturbation sensitivity analysis Of equation (3.3) 


can be performed by computing the total differential of 





equation (5.3) ) that fie, 
oPfa oPfa 
dPfa = [lS a ари rana (3.4) 


Assuming that the SNR is a constant, or equivalently, the 


differential dSNR is zero, dPfa becomes 


ES ا‎ ር De (SET 


OT 


Nees [1+SNR] Pd" Apa (3200 


[ 


where 
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dPfa ~ APfa =  |Pfa -Pfa| 67) 
and 
dpa 2 APd =  |Pd -Pd| (328) 


Dividing both sides of equation (3.6) by equation (3.3) 
yields 

A ۔‎ [1+swr] SES (3.9) 
where APfa/Pfa and APd/Pd represent the fractional error 
between the actual and estimated values of Pía and Pd, 
respectively. 

Such an analysis indicates that the percent error in the 
computed estimate of Pía is linearly related to the percent 
error in the estimate of Pd. The constant of proportionality 
relating the error in the computed estimate of Pía to the 
error in Pd estimate is equal to (1 +SNR) where the SNR is 
taken at the output of the magnitude-square operation. Thus 
for SNR values greater than approximately 1, or 0 dB, a small 
error in the relative frequency estimate of Pd is scaled to 
a larger error in the estimate of Pfa found by using equation 
(3.3). Note that for a 5x5 element planar array anda 
bandwidth of 5 times the fundamental frequency as in the CW 


pulse case, the array element input SNR can be found by 
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equation (2.92), and is roughly 24 dB less than the SNR at 
the output of the magnitude-square operation. That is, a 

0 ав SNR at the output of the magnitude-square operation 
corresponds to an array element input SNR of -24 dB. The 
end result is an estimate of Pfa that diverges wildly from 
the specified Pfa parameter at SNR values over which the 
receiver operates. In addition, equation (3.3) was derived 
using certain assumptions regarding the statistics of the 
input noise that may not be precisely duplicated by the 
pseudorandom noise data generated by the computer program. 
If the assumptions regarding the use of zero mean, uncorre- 
lated, Gaussian noise are not satisfied by the pseudorandom 
noise source, equation (3.3) does not hold, and the algorithm 
would not be suitable for other input noise models. For 
these reasons, the use of equation (3.3) was abandoned in 
favor of an empirically determined fixed number of trials to 
estimate Pd. 

The number of trials needed to estimate Pd was found by 
assuming that the absolute minimum number of trials should 
be the reciprocal of the Pfa parameter. That is, if a Pfa 
of 0.01 is specified, at least 100 trials must be taken to 
allow at least one chance in one hundred of a false alarm 
Occurring. Even though a relative frequency estimate of Pd 
is being computed, Pd is related implicitly to Pfa through 
the decision threshold, and the value of the Pfa parameter 
should be taken into account when attempting to fix the 


number of trials needed to estimate Pd. 
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The number of trials was empirically determined by 
running the receiver simulation program repeatedly with 
increasing multiples of the minimum number of trials, and 
observing the effect on the estimate of Pd. Depending upon 
the value of the Pfa parameter, it was found that four to 
eight times the minimum number of trials would produce 
curves that did not change appreciably as the number of 
trials was increased further. Therefore, the fixed number 
of trials used to estimate Pd was arbitrarily set at 10 x 1768 
ህህ ra of 0.1, and to 5 x l/Pfa for Pfa of 0.01. The smaller 
multiplier for the Pfa of 0.01 became necessary due to limits 
52735 mpucer resources. 

1. Subprogram READY 

The function of subprogram READY is to obtain the 
simulation parameters and the complex envelope output 
electrical signal data from a data file. The data file is 
generated by the ocean communications channel simulation 
computer program developed by Vos and Ziomek [Ref. 3]. That 
is, READY forms the interface between the RCVR simulation 
and the ocean communications channel program. A flowchart 
of READY is shown in Figure 10. 

The simulation parameters are read in first, and 
are stored in the named COMMON blocks: HEADER, SIGNAL, 
ARRAY, and MEDIUM. The named COMMON blocks provide the 
mechanism for communicating simulation parameters into the 


various subprograms. The HEADER data documents the type of 
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Figure 10. Subprogram READS 4 ٤0 


communication channel, and the date the data was generated 
by the ocean communication channel simulation. The SIGNAL 
data includes: the fundamental frequency £5 and the perm 
To; they dumber sor Һелиан Ky си ии ики fo alına 
sample period To; the number of time samples L; the carrier 
frequency Et and the number of zeroes padded to the input 
signal data, if any. The ARRAY data includes: the number 
array elements M and element spacing dy in the x-direction, 
the number of elements N and element spacing dy in the y= 
direction; the direction cosines о. and Vos representing the 


direction of the direct path from the transmit array to the 
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receive array; the direction cosines ор апа Ур. 


the direction the transmit array beam pattern was steered in 


representing 


the communication channel simulation; the depth Ya of the 
center element of the transmit array; and the depth ye of 
the center element of the receive array. The MEDIUM data 
includes the speed of sound Co at the center element of the 
transmit array; the sound-speed-profile gradient g; and 

the speed of sound Cyr at the center element of the receive 
ray. 

The time-sampled, complex envelope, output electri- 
cal signal data is read in next. The data is stored in the 
complex matrix variable YCE vith dimensions L, M and 
277 The maximum values of L, M, and N are limited to 33, 11 
and 11, respectively. 

2. Subprogram SGNLGN 

The function of subprogram SGNLGN is to generate the 
time samples of the complex envelope and to compute the 
energy of the local processing waveform given by equation 
(2.64). A flowchart of Subprogram SGNLGN is shown in Figure 
ii. 

The local processing vaveform 9(2፲.) is synthesized 
from a finite, complex Fourier series with provisions for 
incorporating the estimates of time delay T and doppler 
ehift $ in the total received signal 7757 Accurate esti- 


mates of T and $ are necessary to achieve maximum output 


from the correlator/matched filter detector portion of the 


67 


SGNLGN 


V 













GENERATE LOCAL 
PROCESSING WVFM 
USING FOURIER 
SIGNAL SYNTHESIS 


CONPUTE ENERGY 
IN LOCAL 
PROCESSING 

WAVEFORM 


RETURN 


Figure 11. Subprogram SGNLGN Flowchart 


receiver. That is, the maximum receiver sensitivity is 


obtained when.T = Ta and $ = 


(2.65) and 18 7 


>, as indicated in equations 
For the purpose of this study, the actual doppler 
shift is always set to zero in the transmit signal syn- 
thesizer, and the actual time delay due to the range or 
distance between array 15 set by the system geometry under 
consideration and the reference speed of sound at the 
transmit array. Therefore, the estimate of doppler shift is 
set to zero in Subprogram SGNLGN, and the estimate of time 
delay due to range is computed from the line-of-sight range 


between the center element of the transmit array and the 
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center element of the receive array, and the speed of sound 


Co at izi center eclement” Of the transmit array, or 





(5510) 


where Y او‎ are the coordinates of the center element 


of the transmit array, and (X عو ا‎ ) are the coordinates of 


r 
the center element of the receive array. 

The estimates of T and $ are easily incorporated 
into the local processing waveform of equation (2.64) by 


applying well-known properties of Fourier transforms to 


Lion (2.17) to yield 


. M 3 K A 2 
x(t-1)exp[j2mot] = ) 2-:.:5:55..1 (5711) 
q=-K 


Note that the right hand side of equation (2.64) is just the 
time-sampled form of equation (3.11). 

The complex Fourier coefficients Sy of equation 
(3.11) are identical to those used in equation (2.17) to 
generate the transmit signal in the ocean communication 
channel simulation computer program. Thus, the local 
processing waveform is identical to the transmit signal in 
functional form and total energy content, but is shifted in 
time and frequency by the estimates of range delay and 


doppler shift, respectively. 
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Subprogram SGNLGN also computes the local processing 
waveform signal energy E. for later use in setting the 
decision threshold y. To compute the energy the magnitude- 
square of the complex Fourier coefficients are summed over 
all harmonics aS indicated in equation (2.42) where Cy is 
the q-th harmonic of the local processing waveform. This 
sum is equivalent to the time-average power in the complex 
envelope of the local processing waveform. The energy can 
then be found by multiplying the time-average power by the 
fundamental period To of the local processing waveform. 

3. . SUbprogram AMPWGM 

The function of subprogram AMPWGT is to provide the 
real-valued, amplitude ft +8  ፡ b of the separable 
complex weights S and d, in equation (2.26). A 1607 
of subprogram AMPWGT is shown in Figure 12. 

To generate the rectangular amplitude window, an 
and b. are set equal to 1.0 for all elements (m,n) in the 
receive planar array. A separate subprogram to generate the 
amplitude weights facilitates generation of other forms of 
amplitude windows such as triangular, Hamming, Blackman, 
etc. However, only the rectangular window is used in this 
study. 

4. Subprogram PHSWGT 

The function of subprogram PHSWGT is to generate the 

E factors Om and on of the separable complex weights 


Cm and a, in equation (2.26). A flowchart of subprogram 
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Figure 12. Subprogram AMPWGT Flowchart 


PHSWGT is shown in Figures 13a and l3b.  PSHWGT can be 
programmed to compute phase factors that compensate, or 
remove the effects of, system geometry and deterministic 
medium wave propagation effects.  PHSWGT can also introduce 
random noise in the phase factors for study purposes. 

The phase corrections for system geometry are com- 
puted using equations (2.27) and (2.28), where the direction 
cosines u and v are selected such that receive array beam 
pattern is aimed at the transmit array along the direct path 
from the receive array to the transmit array. That is, if 


^ 


Hin equation (2.5) represents the direction of the direct 


path from the transmit array to the receive array, u and v 


71 


۲۸ 


SET LOOP COUNTER 
LIMIT VARIABLES 
TO INITIAL VALUES 









SET DIRECTION 

COSINES U,V TO 
DESIRED DIRECTION 
U=UB,UR V=VB,VR 















COMPUTE SPATIAL 
FREQUENCIES FOR 
X AND Y 
COORDINATES 








SET GEOMETRIC 

PHASE CORRECTION 
FACTORS TO 

ZERO 







COMPUTE GEOMETRIC 
PHASE CORRECTION 
FACTORS 


Figure 13a. Subprogram PHSWGT Flowchart 
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Figure 13b. 






are set equal to “Uy and “Vo So that the receive array 
beam pattern points in the . direction. 

The deterministic medium phase correction factor 
is found by using equation (2.32) to compute the phase shift 
of the signai due to propagation through an inhomogeneous 
medium, and then negating the result. The random medium 
effect can also be computed from a closed form expression, 
but its use in generating phase correction factors is not 
the subject of this study, and will not be discussed. The 
total phase correction factor for the y-direction is 
obtained by adding the system geometry and deterministic 
correction factors: 

5. Subprogram ARYPRO 

The array processor subprogram ARYPRO uses DFT and 

IDFT algorithms to: 


- generate the spectrum of the input electrical signal 
data at each element, 


- correct the phase of the spectral components to co- 
phase the signals at all array elements, and 


- inverse transform the signal spectrum at each array 
element to recover the cophased signal data. 


A block diagram of the subprogram ARYPRO is shown in Figures 
14a and 14b. 

The output of the array processor is a time-sampled, 
complex envelope signal representing the sum over all array 
elements of the signals at each element. The defining 
equation for the DFT is used instead of a fast Fourier 


transform algorithm because the maximum number of time 
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Figure 14b. Subprogram ARYPRO Flowchart 


samples permitted (33) is insufficient to achieve a measurable 
improvement in speed of execution [Ref. 6:pp. 151-152]. 
Additionally, use of the DFT equation relates the structure 
of the program directly to previous work by Ziomek [Ref. 2] 
and Vos [Ref. 3]. 
6. Subprogram AWGN 

The function of subprogram AWGN is to generate one 
sample of an uncorrelated, Gaussian process with arbitrary 
mean and variance. A flowchart of AWGN is shown in Figure 15. 

AWGN is based on the International Mathematical 
Subroutine Library (IMSL) FORTRAN pseudorandom number 
generator routine GGNQF.  GGNQF is a function subprogram 
that returns One zero mean, unit variance, Gaussian, or 
N(0,1), pseudorandom number with each call to the function 
Subprogram. The zero mean, unit variance pseudorandom 
number X is then scaled with the desired mean H. and the 


desired standard deviation om using the relation 


Xf = 9 7. (25-12) 


vhere the standard deviation is computed from the variance 


as 


o = O (SES) 


The desired mean and standard deviation are passed 


to AWGN as arguments from the main program. The mean is 


TI. 


DRAW ONE NUMBER 

FROM NORMAL(0,1) 

RANDOM NUMBER 
GENERATOR 


SCALE NUMBER 

VİTH DESİRED 

MEAN AND STD, 
DEVIATION 





Figure 15. Subprogram AWGN Flowchart 


always set to zero for the purposes of this study and the 
Standard deviation is found from equation (3.13) where the 
variance represents the desired power level in either the I 
or Q channel noise signal. 

The function subprogram GGNQF internally recomputes 
a new seed value for subsequent calls from a seed provided 
on the first call to the function. Since the first seed 
value is set aS a parameter in the top level program, and 
is passed as an argument through AWGN to GGNOF, the pseudo- 


random sequence always follows the same pattern each time 
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the receiver simulation is run. Using the same pseudo- 
random sequence for each simulation run is essential if 
comparisons between different runs are to be made from the 
receiver simulation output plots. 
7. Subprogram INTGRT 

The function of subprogram INTGRT is to compute, 
HNcCegquation (2.68), the approximation of the correlation 
Nu ral grven in equation (2.67). A flowchart of subprogram 


NEP 15 Shown in Figure 16. The complex product of the 


INTGRT 


SUM ALL TERMS 

OF COMPLEX- 
VALUED INPUT 
SEQUENCE 















MULTIPLY SUM 
BY TIME 
INCREMENT 

Ts 


RETURN 


Figure 16. Subprogram INTGRT Flowchart 





total signal and noise, and the local processing waveform 


is computed in the top level program. The resulting 


zə 


complex-valued Sequence is passeli to INTCR иһаси сәтири. 
the correlator our ue Ë by using a trapezoidal integration 
algorithm to approximate the integral of equation (2.67). 
A separate subprogram allows easy implementation of other 
numerical integration algorithms, if desired. 
8. Subprogram WRITER 

The function of subprogram WRITBL is to simply 222 
vide tabular output of the simulation parameters and selected 
information generated during the execution of the program. 


A flowchart of subprogram WRITBL is shown in Figure 17. 


WRITBL 


WRITE 
SIMULATION 
PARAMETERS 

TO OUTPUT FILE 









Figure 17.  Subprogram ۷۸ Flowchart 


Besides the data read into the COMMON storage blocks 


by subprogram READY, the information output to the data file 
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in table form includes the array element input SNR at which 
the Pd is estimated, the estimate of Pd, an estimate of Pfa 
obtained using equation (2.82), the SNR at the output of the 
array processor, the array gain and the SNR at the output of 
the magnitude-square operation. The tabular output is 
intended primarily to provide a convenient means of testing 
a ram modifications. It is not intended to be the primary 
simulation output, and will not be discussed further. 

7 Subprogram PDPLOT 

The function of subprogram PDPLOT is to convert the 
numeric data generated by the receiver simulation into 
graphic form. A flowchart of the subprogram PDPLOT is shown 
in Figure 18. The plots generated from the numeric data 
ExcEobtarned through the use of standard DISSPLA graphics 
library subroutines. 

The primary output of the receiver simulation is a 
plot of the estimate of Pd versus the array element input 
SNR with a fixed value of Pfa as a parameter.  PDPLOT also 
computes a theoretical value of Pd using the equation (2.82) 
where the SNR at the output of the magnitude-square operation 
is related to the array element input SNR through equation 
E91). 

The curve representing the estimated Pd is plotted 
through the estimated Pd data using the DISSPLA least squares, 
emic spline, curve fitting plot routine SMOOTH. The value 


of Pd obtained from equation (2.82) is plotted using the 
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Figure 18.  Subprogram PDPLOT Flowchart 
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DISSPLA cubic spline, interpolating polynomial routine 
SPLINE, and appears as a dashed line on all plots generated 
by the receiver simulation. That is, the curve for the 
calculated value of Pd is an interpolating polynomial that 
passes through the calculated Pd data points, while the 
curve for the estimated value of Pd is allowed to pass 
within some arbitrary, small offset of the estimated Pd 
data point. For the purpose of this study, the maximum 
offset allowed was 10 percent of the data point value. 

This value of offset provided a reasonable compromise be- 
tween obtaining a relatively smooth fit of the plotted curve 
without diverging excessively from the data points. Both 
the estimated and calculated Pd data points are plotted at 


the same array element input SNR. 


277 MODEL VERIFICATION 

Nerlification of the receiver model involved two tasks, 
characterization of the pseudorandom number generator noise 
source, and a comparison of the data generated by the re- 
ceiver simulation program with the results predicted by 
the theoretical development discussed in Section II. Why 
simulate a test case for which a theoretical model exists? 
The main reason is to gain some confidence in the results 
generated by the simulation when the inputs are such that a 
theoretical model does not exist or is mathematically 


iüractable. 
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41. . ሮከ፲ውር!ሮት1ሪክ ሰሰጠ፡።፡ዚሀ ፡፡ገርጆር ፡ርየ e 
The capability to generate noise samples consistent 

with the noise description used in developing the receiver 
equations is central to verification and usefulness of the 
computer simulation. Thus, the first step is to ascertain 
the statistical properties of the noise source used in the 
computer program. Recall that the assumption required in 
developing therelations between Pd, Pfa, decision threshold 
and SNR of the magnitude-square of the correlator output 
involved the use of zero mean, Gaussian noise which is uncorre- 
lated in both spatial and temporal coordinates. The noise 
Source was tested by taking a sequence of noise samples for 
both I and Q channels in manner identical to that used in 
the receiver simulation. That is, the subprogram AWGN 
was embedded in a test program to ascertain the statistical 
properties of the complex envelope of the noise signal. The 
following tests were applied to the I and Q channel sample 
sequences to verify agreement with the noise process 
assumptions: 

- sample mean 

- sample variance 

- histogram 

- autocovariance, and 

- estimated power spectral density. 
The test algorithms were obtained from standard IMSL 


procedures. 
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Figures 19 and 20 show the results of the sample 
mean, variance and histogram calculations for the noise 
processes in the I and Q channels, respectively. The sample 
mean and variance calculations approximate the zero mean, 
unit variance assumptions quite well. The histogram 
calculations indicate the distribution of the samples is 
roughly Gaussian. The smooth curve plotted through the 
histogram represents the exact Gaussian distribution for the 
size of the sample window used. The greatest deviation from 
Gaussian appears near the mean value which unfortunately is 
where most of the sample values lie. Thus, some difficulty 
in achieving a perfect correlation between simulated per- 
formance and theoretical results was anticipated, and in 
fact some deviation from theory at large values of Pd did 
Soccur. 

The space-time correlation properties of the sample 
sequence was measured by computing the autocovariance func- 
tion of 2000 samples over a total time (or space) displace- 
ment (lag) of 1000 samples. The autocovariance was computed 
using the IMSL routine FTAUTO. The correlation of adjacent 
samples, whether one or two or 1000 samples apart, was 
found to be remarkably small. The autocovariance function 
of the I and Q channel noise is shown in Figures 21 and 22, 
respectively. Because of the manner in which the samples 
are drawn in the simulation, the distinction of a sample 


being assigned to a particular channel at a particular 
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element at a particular time is immaterial. The auto- 
covariance of the samples indicates that noise samples can 
be generated which are independent of sample time, array 
element location, and receiver channel. 

The estimated PSD was computed primarily to rein- 
force the results obtained from the autocovariance compu- 
tation. The IMSL routine, FTFPS, used to compute the 
PSD of the sample sequence implements an algorithm similar 
to that due to Welch [Ref. 7:pp. 553-554] in which the term 
sample record iS partitioned into contiguous subrecords. 
Each subrecord is иә weighted with a triangule window 
function. Then, a periodigram of each amplitude weighted 
sample subrecord is computed. The resulting subrecord 
periodigrams are averaged over all subrecords in an effort 
to reduce the variance of the estimate of the PSD. The 
estimate of PSD for the I and Q channel noise is shown in 
Figures 23 and 24, respectively. The estimated PSD is 
approximately flat for both the I and Q channels over the 
sampled frequency range of 0 to m radians, or one sample 
frequency period. This tends to support the results obtained 
from the autocovariance computation that the noise samples 
are indeed uncorrelated. 

2. Verification of the Output Data 

The primary output data from the receiver simulation 

is a plot of an estimate of Pd for various values of input 


SNR at a given Pfa. As shown in equation (2.82), a closed 
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form relation exists between the Pd, Pfa and correlator 
output SNR. Furthermore, the correlator output SNR can be 
related to the array input SNR when the signals at all array 
elements are precisely cophased by equation (2.91). 
Therefore, a plot of the theoretical performance computed 
using equations (2.82) and (2.91) can be superimposed on 
the plot of the receiver simulated performance. The curve 
representing theoretical performance provides verification 
of the simulated performance when input signal parameters 
and phase weighting allow precise signal cophasing to 
occur. The computed performance curve also provides a 
baseline to evaluate simulated performance when the input 
noise or receiver operating conditions differ from the 
BEDNNUptrons used in formulating the receiver model. In 
all the plots, the theoretical performance is shown as a 


dashed curve. 


ds 
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This section presents the results of three receiver 
simulation case studies and analyzes the simulation output 
in each case. The three cases considered are: 

- Case TEST, in which the output electrical signal data 
at each array element is produced by a test signal 
generator computer program, and is free of communi- 
cation channel effects such as attenuation, phase 
shifts due to wave front refraction, time delay due to 
range separation, and the transmit array beam pattern. 

- Case HMG1, in which the output electrical signal data 
at each array element is produced by the ocean 
communication channel simulation computer program, 
and includes the effects of path attenuation, phase 
shifts due to system geometry, time delay due to range 
separation, and the transmit array beam pattern. 

- Case INHMG1, in which the output electrical signal 
data at each array element is produced by the ocean 
communication channel simulation computer program, 
and includes not only the effects given in case HMGl, 
but also contains phase shifts due to refraction of 
the wavefront along the propagation path. 

The results of these Case studies will be analyzed by 
providing a brief description of the transmit signal 
waveforms used to generate the array element output elec- 
trical signal data, by listing the system parameters that 
distinguish the test cases, and finally, by interpreting 
the plots of Pd versus array element input SNR for each 


case study. 


A. TRANSMIT WAVEFORMS 
The analysis and theoretical development of the receiver 


simulation in Section II makes no assumption regarding the 
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functional form of the transmit signal. The receiver 

should perform equally well regardless of the type of trans- 
mit signal provided that the correlator/matched filter uses 
a replica of the transmit signal as the local processing 
waveform. To verify this hypothesis, and thereby test the 
integrity of the receiver simulation and the validity of the 
assumptions used in generating the model, more than one type 
of transmit signal was used to test receiver simulation 
performance. The use of multiple transmit signal waveforms 
also provides a broader base from which to draw conclusions 
regarding the concept of model-based signal processing. For 
the purpose of this research, two types of transmit wave- 
forms were used in producing the Pd versus array element 
input SNR plots in each of the three case studies. 

A rectangular envelope CW pulse and a rectangular 
envelope LFM pulse were selected as transmit signals. These 
two particular forms were chosen for several reasons. First, 
the CW pulse and LFM pulse waveforms were considered to be 
typical transmit signals used in SONAR and acoustic signal- 
ing systems. Second, the CW pulse provides a simple case 
of amplitude modulation while the LFM pulse provides an 
example of an amplitude and angle modulated waveform with 
considerably different spectral characteristics. Finally, 
either signal can be readily synthesized from a finite, 


complex, Fourier series whose coefficients Em may be calculated 
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from a closed form expression for the frequency spectrum 
of the time domain pulse characteristics. 

The rectangular envelope CW pulse is discussed first. 
The pulse repetition frequency is the same as the funda- 
mental frequency . of the finite frequency spectrum from 
vhich the pulse is synthesized. The pulse duty cycle vas 
arbitrarily taken to be 0.5 yielding a pulse vidth of Т./2 
seconds where To = 1/Е2. The complex, Fourier series 
coefficients used to synthesize the complex envelope of the 
CW pulse are obtained from a closed form expression for the 
complex-valued continuous spectrum. To obtain the Pour ii 
coefficients, the closed form expression for the continuous 
spectrum is evaluated at the discrete frequencies Ед, and 
the resulting complex value is divided by the fundamental 
pulse period T5: The index q takes on integer values 
between -K to K where K is the maximum number of harmonics 
in the finite Fourier series used to synthesize the CW 
pulse. The continuous spectrum of the CW pulse is the 
familiar sin(x)/x form obtained by taking the Fourier trans- 
form of the rectangular pulse shape where the pulse width 
is one-half the pulse period. The following specific 


transmit signal parameters were used in all CW pulse 


simulations: 
- Amplitude, A = 40.0 
= Duty Cycle; D 2 95 


200 0 


Fh 
| 


- Fundamental Frequency, 
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- Number of harmonics, KMAX = 5 


- Harmonic values. Co = 2000000 ëexp[50°] 
C = CJ = 12.33240 exp | 70] 
Cc 2” c, = 0.000000 exp[j0°] 
وره‎ Z C} = 4.244132 exp[j31809?] 
TES 0.۰00000 sex D | | 
CE e. — 2.546479 exp[j0?] 


The LFM pulse, or linear frequency modulated pulsed 
carrier waveform, is discussed next. The complex Fourier 
coefficients used to synthesize the LFM pulse are found 
using a procedure similar to that used for the CW pulse 
except the closed form expression for the complex-valued 
continuous spectrum of the LFM pulse was found by using the 
method of stationary phase. Officer [Ref. 8:pp. 67-68] 
describes the method of stationary phase as does Papoulis 
[Ref. 9:pp. 267-273] who also provides a complete descrip- 
tion of the LFM waveform. The following transmit signal 


parameters are used in all LFM pulse simulations: 


- Amplitude, A = 40.0 

x Duty Cycle, D - 0.8 

- Phase Deviation Constant B = 2356.2 radians/volt 

- Fundamental Frequency, f = IO Hz 

- Number of harmonics, KMAX = 3 

- Harmonic values, Co = 14.60593 exp[j45°] 
С̧ј = сј = 14.60593 ехр11212) 
Cc” Q 14.60593 exp13309”) 
С̧ј = с̧у = 14.60593 ехр131892| 


It should be emphasized that these complex, Fourier 
series coefficients which are used in the ocean communica- 
tion channel simulation computer program to synthesize the 
complex envelope CW and LFM pulse transmit waveforms are 
also used in subprogram SGNLGN to produce the local processing 
waveform for the correlator/matched filter detector section 


of the receiver simulation computer program RCVR. 


be +0 PE ST 

The input signal data for case TEST is obtained from a 
Separate computer program written to synthesize array 
element output electrical test signals from a finite Fourier 
Series expansion of the test signal waveform. The technique 
is similar to that used in generating the receiver simula- 
tion local processing waveform in subprogram SGNLGN.  How- 
ever, appropriate phase shifts are added to the Fourier 
coefficients to vary the direction of arrival of the plane 
wave incident on the receive array, and to incorporate the 
time delays due to spatial separation, d. and dr of the 
array elements. Furthermore, the time delay due to range 
and the doppler shift of the test signals are both set to 
zero to prevent signal loss in the correlator/matched filter 
due to a mismatch between the total array output signal and 
the local processing waveform. In this way, array element 
output electrical signals are obtained that exclude effects 
due to propagation of the sound wave through the ocean 


medium, or spatial modulation of the signals due to the 
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transmit array beam pattern. Since the medium effects are 
excluded from the case TEST array element signals, the 
carrier frequency Lo is set to zero, and the array element 
Spacing is increased to keep the interelement spacing d. 
and dy equal to one-half the shortest wavelength in the 
transmit signal. Note that this produces different inter- 
element spacings for the CW and LFM pulse situations. 

Signals generated in this fashion allow the receiver 
output to be analyzed separately from the communication 
channel. Essentially, it is equivalent to having a signal 
generator that provides a test signal to measure receiver 
performance. The results obtained from case TEST, when 
combined with the plot of theoretically predicted performance, 
provides a baseline with which to analyze the receiver 
7 when the ocean communication channel simulation data 
1s processed. 

The system parameters defining case TEST are as follows: 

— Array Parameters 
Number of array elements, MN 25 


0.7375 meters 


ر0 
| 

О, 
1 


Array Element Spacing, CW 


О, 
|| 
О, 
| 


LFM 24.58 meters 


- Medium Parameters 


Speed of Sound, əm 1475 meters/second 
Actual time delay Tes 0.0 seconds 
Actual doppler shift Фд = 050 Hertz 


Four plots representing receiver performance were 


generated for case TEST. Phase weighting is done to 


99 


compensate for the direction of arrival of an acoustic plane 
wave incident upon the array elements. The direction of 
arrival can be set to any arbitrary direction by the test 
Signal generator computer program. However, the direction 
was chosen to be the same as that due to the system geometry 
used in the case studies HMGl and INHMGl. Phase weighting 
for geometry is indicated by the state of the logical 
variable STEER. The estimates of time delay due to range 
and doppler shift were both set to zero in subprogram SGNLGN 
to maximize the output of the correlator/matched filter 
detector. Two of the plots, Figures 25 and 26, show receiver 
performance for a rectangular envelope CW pulse. Figure 25 
was produced using a Pfa of 0.1, and Figure 26 shows the 
effect of a Pfa of 0.01. The remaining two plots generated 
for case TEST, — 27 and 28, show receiver performance 
for a rectangular envelope LFM pulse at a Pfa of 0.1 and 
0.0l, respectively. Since the input electrical signal data 
does not incorporate channel effects, phase weighting for 
deterministic and random medium effects was not done for 
case TEST. 

Several observations regarding model performance can 
be made by analyzing the plots. Generally for all plots, the 
measured performance of the simulation (solid curve) follows 
the predicted theoretical performance (dashed curve), but 
the agreement is not exact. That is, the estimate of Pd 


follows the same trend as the predicted Pd, but the 
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K I Or from a low Pd to a Pd of 1 is more sensitive to 
changes in the array element input SNR. A comparison 
between the plots shows the Pd to be a function of the Pfa 
as predicted by equation (2.82), and for small values of 
SNR, the estimated Pd approaches the Pfa parameter used by 
the simulation in generating the plots. This effect is 
also predicted by equation (2.82) since at SNR values much 
less than one, the input signal is essentially off, and only 
the noise is producing signal detections. A comparison 
between the plots for the CW pulse and LFM pulse does show 
some minor difference in receiver performance, but the same 
general trends seem to hold regardless of the functional 
form of the transmit signal. - 

The differences in the measured performance between the 
CW pulse and the LFM pulse simulations can be attributed in 
part to the difference in the maximum number of harmonics 
used to synthesize the transmit signals. As shown by 
equation (2.22), the number of time samples L is related to 
the maximum number of harmonics K used to synthesize the 
transmit signal. The difference in the number of time 
samples causes different subsequences of noise samples to be 
drawn from the noise generator subprogram AWGN. The 
different noise sample subsequences will cause small varia- 
tions in the relative frequency estimate of Pd. Thus, an 
exact comparison of the noise performance cannot be made 


across different types of transmit signals when the number 
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of harmonics used to create the transmit signal are not the 
same. However, the real significance of this effect is that 
even with different noise sample subsequences in the I and 

Q receiver channels at the array elements, the measured 
performance conforms closely to the performance predicted by 
theory under the very ideal assumptions of precise signal 
cophasing, and zero-mean, white, Gaussian input noise. 

For these reasons, it can be inferred that the basic 
receiver model is valid, but some difficulty may exist in 
generating precisely white, zero-mean, Gaussian noise within 
the simulation program. Indeed as indicated in Section 
III.B.1, the noise source used in the receiver simulation 
does not generate precisely zero-mean, white, Gaussian 
noise. Furthermore, the measured Pd is a relative frequency 
estimate of the actual receiver performance, and may also 
contribute some error in the simulation results. 

Considering these approximations to the assumptions used 
in the receiver model development of Section II, exact 
agreement should not be expected even in the case where the 
communication channel effects are excluded. However if 
only relative changes or effects are to be observed, moderate 
disagreement between the theory and simulated receiver 
performance can be negated by comparing only differences in 
the estimate of Pd for each of the case studies. A com- 
parison of the differences in the measured receiver perfor- 


mance will indicate if, and to what degree, the effects due 
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to the physics of acoustic wave propagation can be compen- 
sated for by model-based signal processing. The results of 
case TEST then can be viewed as a validation of the receiver 
model and the applied programming used to implement the 


model. 


C. CASE HMGlI 

The input signal data for case HMGİl and case INHMGl 15 
produced by the ocean communication channel simulation 
program due to Vos and Ziomek [Ref. 3]. For case HMGl, the 
OCean is considered to be a homogeneous medium with respect 
to the speed of sound propagation through the water. That 
is, the speed of sound is identical at the transmit and 
receive arrays, and at all points in between. In this model 
of the ocean medium, refraction or ray bending of the sound 
wave does not occur. Thus, only phase weighting for geometry 
is required, and only those results will be presented. The 
homogeneous case provides a needed baseline with which to 
judge the effects of model-based signal processing: when 
inhomogeneous case data is studied. 

Since case HMGİl includes the effects of the ocean com- 
munication channel, a carrier frequency is used to convey 
the modulation through the channel. In case HMGl and case 
INHMG1 the carrier frequency was set to 5.0 kHz. The array 
interelement spacing was adjusted accordingly to maintain 


the one-half wavelength spacing needed to prevent grating 


ın. 


lobes in the array beam patterns. Note that the interelement 
Spacing was set using the highest frequency expected for 
either the CW pulse or the LFM pulse waveforms, and was 
kept constant for both the HMGl and INHMGIl case studies. 
The highest frequency occurred for the CW pulse case, and was 
equal to 6000 Hz. The estimates of the time delay and 
Doppler shift in subprogram SGNLGN were set equal to the 
actual values to maximize the correlator output. 

The system parameters defining case HMGl are as follows: 


- Array Parameters 


Number of array elements, ከ1 = ዝ = 5 
Array Element Spacing, dy š: . - 0.1229 meters 
- Medium Parameters 
Speed of Sound, cn 1475 meters/second 
Actual line of sight 
time delay Ta = 2.033898 seconds 
Actual Doppler shift On = 0.0 Hertz 
- System Geometry (See Figure 1) 
Depth of Transmit Array y 1000.0 meters. 
Depth of Receive Array y 2500.0 meters 
Cross Range CIE C MEE 500.0 meters 
Line of sight range |r-r.| - 3000.0 meters 


Eight plots were generated for case HMGl. The first four 
plots, Figures 29 through 32, were produced using the 
rectangular envelope CW pulse transmit signal. The remain- 
ing four plots, Figures 33 through 36 were produced by 


the rectangular envelope LFM pulse transmit signal. Figures 
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29, 31, 33 and 35 show receiver performance when no phase 
weighting is done. These plots indicate that performance 
will be very poor when the output signals from the array 
٣٣۰۰٠۰٣۶۰ are not cophased. Again, it is significant that 
the Pd in this case is approximately equal to the value of 
the Pfa parameter used in obtaining the plot. That is, 
the total array output signal must be very small because 
Of destructive interference of the array element signals, 
and the only detections appear to be due to the noise. 

Figures 30, 32, 34 and 36 show performance when phase 
weighting is done to compensate for system geometry. These 
plots indicate receiver performance can be made to approxi- 
mate the theoretical predicted performance if the individual 
array element signals can be cophased to provide constructive 
interference at the output of the array processor. 

Figures 29, 30, 33 and 34 show receiver performance at 
0 11١33 or 0.1. Figures 31, 32, 35 and 36 show the effect of 
decreasing the Pfa to 0.01. The effect of decreasing the 
Pfa on the simulated performance follows the theoretically 
predicted performance and provides further confirmation that 
the computer implementation of the model is accurate. 

The important conclusion from case HMGl is dramatic 
277 venent in receiver performance is possible provided 
that enough information about system geometry exists to 
generate phase weights that will properly cophase the array 


element output signals. 


18800 


İb: CASE INHMGI 

The input signal data for case INHMGl is produced in 
the 257 manner as case HMGl. The transmit signal vaveforms, 
system geometry, and transmit and receive arrays are identi- 
cal. The only difference betveen the tvo cases is the 
model of the ocean medium. For case INHMGl, the ocean is 
considered to be an inhomogeneous medium vith respect to the 
speed of sound propagation through the vater. The speed of 
sound is assumed to vary linearly vith depth (or y-coordinate) 
only. That is, the speed of sound is not identical at the 
transmit and receive arrays. In this space-variant model 
of the ocean medium, refraction or ray bending of the sound 
wave does occur, and the model of the ocean medium due to 
Ziomek [Ref. 2] predicts what the phase shifts in array ele- 
ment output electrical signals will be. Using this knowledge 
of the physics of acoustic wave propagation, phase weights 
are generated that attempt to compensate for the refraction 
of the acoustic wave allong che ray path” and hoe ии 
geometry. To determine the effect of this model-based signal 
processing algorithm, the inhomogeneous data was also 
processed using phase weights that compensate just for sys- 
tem geometry. The difference in simulated receiver perfor- 
mance between these — situations vill graphically shov the 
effectiveness of the model-based signal processing approach. 

The system parameters defining case INHMGİ are as 


follovs: 


18 


۰۰۰۲۴۹۲ ۲۹٣۹۶۹۲۲۷٠٢5 


Number of array elements, М = Ы = 5 
Array Element Spacing, a, = ду کے‎ 0.1229 meters 
- Medium Parameters 
Speed of Sound, c s 1475 meters/second 
Gradient g = .017 
Actual line of sight 
time delay Uu 2.033898 seconds 
Actual Doppler shift Фд = 0 0SHertz 
- System Geometry (See Figure 1) 
Depth of Transmit Array "ES 1000.0 meters 
Depth of Receiver Array ےر‎ 2500.0 meters 
Cross Range 20556” 500.0 meters 
Line of sight range in = 3000.0 meters 


Twelve plots were generated for case INHMG1. The first 
six plots, Figures 37 through 42 show receiver performance 
for a rectangular envelope CW pulse. The remaining six 
plots, Figures 43 through 48, are the result of a rectangu- 
lar envelope LFM pulse transmit Signal. Differences in 
receiver performance due to the form of the transmit signal 
are meaSurable, but not significant. 

Figures 37, 40, 43 and 46 were generated when no phase 
weighting is applied to the array output electrical signals. 
Again, these figures indicate that the array element output 
electrical signals must be cophased to provide any useful 


detection capability in the receiver. 


JESUS, 
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Figures 38, 41, 44 and 47 show the receiver performance 
when oniy the system geometry is compensated for by the 
phase weights used in the array processor. These figures 
indicate that in the presence of the inhomogeneous ocean 
medium, phase weighting to compensate for system geometry 
will not achieve the receiver performance predicted by 
theory. That is, traditional beam steering is not suffi- 
cient for maximum receiver performance, and a large margin 
for improved performance exists when the effects of the 
Ocean medium on the received signal can be predicted. Note 
also that the degradation in performance is consistent for 
the different transmit waveforms and values of Pfa. 

Figures 39, 42, 45 and 48 show the effect on receiver 
performance when phase weights are computed based on both 
the system geometry and the refraction of the acoustic wave 
in the ocean medium. The state of the logical variable 
DMEDIA indicates if the phase weights contain the correction 
for deterministic, inhomogeneous medium wave propagation 
effects. As expected, the simulated receiver performance 
conforms to the performance predicted by theory, and is 
nearly identical to that of case TEST when all medium 
effects have been corrected for. 

However, the really significant result is seen by com- 
paring Figures 38, 41, 44 and 47 with Figures 39, 42, 45 and 
48. The only difference between these plots is in use of 


phase weights that compensate for the phase shifts due to 


ے3ل 


propagation through an inhomogeneous ocean medium, and the 
model of the inhomogeneous ocean medium is based on the 
physics of acoustic wave propagation. If one assumes a 
marginal detection probability of 50% (Pd = 0.5), the 
improvement in receiver performance due to the application 
of a model-based signal processing approach for this particu- 
lar simulation geometry and ocean medium characteristic 
varies from 14 dB to 18 dB depending on the type of 

transmit waveform or value of Pfa. The conclusion is that 
the physics of wave propagation can be used to develop 
Signal processing algorithms having significant impact on 
the performance of a receiver designed to process complex 
envelope, electrical signals from an array of point source 


elements. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The original objectives of this research have been accom- 
plished. The computer simulation of a correlator/matched 
filter receiver to process data from an array of መ ሌም 
acoustic transducers has been developed and validated. The 
concept of model-based signal processing was applied to the 
development of an array signal processing algorithm, and was 
shown to have a marked impact on the capability of the 
receiver to detect the presence of a signal in zero-mean, 
additive, white Gaussian noise. 

The computer Simulation of the receiver was validated by 
a direct comparison between the measured performance of the 
Simulation and the performance predicted by theory when all 
array element output electrical Signals are precisely 
cophased. The agreement between predicted and measured 
performance was close but not exact Since the noise assump- 
tions made in developing the theoretical performance could 
not be precisely duplicated by the computer simulation. 
However, the close agreement with theoretical performance 
was considered to be a validation of the receiver simulation 
667101011067 program; 

Test cases were designed and run that permitted relative 
changes in receiver performance to be measured. In this 
way, the effectiveness of the model-based signal processing 


approach could be quantified. For the simulation parameters 
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and transmit signals used in this study, compensating for 
5 refractıon of the acoustic field in the ocean medium 
was found to improve receiver performance over traditional 
beam steering by at least 14 dB when measured at the point 
of 50$ probability of detection. 

Suggestions for further study of the receiver model 
include: 


- Obtain, or develop, a pseudorandom noise generator 
that more closely approximates the noise model 
assumptions implicit in the derivation of the 
theoretical receiver performance, and revalidate 
the computer simulation. 


- Simulate more realistic noise processes, such as 
colored noise, and determine the impact on receiver 
performance. 


- Try decision rules other than Neyman-Pearson, such 
as the minimum average probability of error criterion 
to measure probability of error performance for an 
underwater, acoustic data communication system. 


- Experiment with other functional forms for the 
transmit signal, and measure the effect on perfor- 
mance when errors in the estimate of time delay and 
Doppler shift exist. 
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